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Abstract 



A random distribution of inert sites is introduced in the Ziff-Gulari-Barshad 



model to study the phase transitions between active and poisoned states. The 
adsorption of CO and O2 molecules is not possible at the position of the inert 

a 

sites. This model is investigated in the site and pair approximations, as well 

-i 

as through Monte Carlo simulations. We determine the mean coverages of the 

o . 

elements as a function of the dilution and show that the continuous transition 
between the active and O-poisoned state is slightly affected by moderate val- 

> 

£2 ■ ues °f dilution in the pair approximation and in the simulations. On the other 

OO ■ hand, from the analysis of the hysteresis curves, the transition between the 

O, 

active and CO-poisoned states changes from first-order to continuous as one 

increases the concentration of inactive sites. The observed transition in the 

site and pair approximations is always of first-order nature. We also found 

C . the lines of transition and spinodal points as a function of the concentration 

of inert sites. Finally, the production rate of CO2 is calculated as a function 
r> ■ 

^ , of the dilution of sites. 

b . 

PACS number(s): 05.70.Ln, 05.70.Fh, 82.65.Jv, 82.20.Mj 
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I. INTRODUCTION 



The study of nonequilibrium phase transitions is a topic of growing interest due to its 
application to a variety of complex systems 1 ' 2 : contact process, domain growth, catalysis, 
phase separation and transport phenomena. Although there is no general theory to account 
for nonequilibrium model systems, in recent years some progress has been achieved in under- 
standing the stationary states of these systems employing approximate analytical methods 
and simulations. Some rigorous mathematical questions concerning the phase transitions of 
these complex interacting particle systems can be appreciated in the books of Liggett 3 and 
Konno 4 . 

In this paper we focus our attention on the phase transitions observed in the surface 
reaction model proposed by Ziff, Gulari and Barshad 5 (ZGB), which describes some kinetic 
aspects of the oxidation of CO over a catalytic surface. In particular, here we consider a 
modified version of the ZGB model, where we incorporate a random distribution of inert 
sites on the catalytic surface. The original ZGB model is an irreversible lattice model for 
surface reactions based on the Langmuir-Hinshelwood mechanism, where the reactants must 
be adsorbed before reacting. The steps used to describe the ZGB model ( a lattice Markov 
process) are the following: Molecules of CO and 2 from a gaseous phase can be adsorbed 
onto the sites of a regular square lattice of identical sites. These molecules arrive at the 
surface according to their partial pressures in the gas mixture, that is, the probability of a 
CO molecule arriving is y co and (1 — y co ) for the 2 molecule. The CO molecule requires only 
a single vacant site to be adsorbed, while the 2 is adsorbed if it finds a nearest-neighbor pair 
of empty sites. Upon adsorption, the 2 molecule dissociates and the two free O atoms can 
react independently. If, after an adsorption step, a nearest-neighbor CO — O pair appears 
on the lattice, they immediately react, forming a C0 2 molecule that goes to the gas phase, 
leaving two empty sites on the lattice. Therefore, in this adsorption controlled limit, only a 
single parameter (y co ) is sufficient to describe the dynamics of the model. 

The simulations performed by Ziff and co-workers have shown that the system exhibits 
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two phase transitions between active and poisoned states: for y co < yi, an O-poisoned state 
is formed, while for y co > y 2 the lattice is poisoned by CO. For yi < y co < y 2 a reactive 
steady-state is found, in which a nonzero number of vacant sites is present in the lattice. At 
yi the transition is continuous, whereas at y 2 the transition is of the first-order type. Using a 
mean field theory, Dickman 6 qualitatively reproduced the phase diagram of the ZGB model 
and showed that, at the level of site approximation, only the first-order transition appears. 
However, employing the pair approximation, both continuous and first-order transitions are 
obtained. 

We are interested on the effects of inert sites on the phase transitions of the ZGB model. 
We have investigated in detail the dependence of the phase transitions on the concentration 
of inert sites. This problem presents some experimental interest in the automobile industry, 
where lead particles are deposited over the catalyst during the exhaust of the gases after 
combustion. This affects the efficiency of the catalytic surface due to the pinning of these 
lead particles on the surface, forbidding the adsorption of CO and 2 molecules at the lead 
positions and reducing the reaction paths. Hovi and co-workers 7 , have studied by computer 
simulations the effect of preadsorbed poison and promoters on the irreversible ZGB model. 
They calculated the coverage of species as a function of the concentration of inert sites for 
a wide range of values, finding the interesting result that the first-order transition changes 
to a continuous one at a critical value of the concentration. Cortes and Valencia 8 have also 
reported some results concerning random impurities distributed over the catalyst, in which 
they observed the change of the first-order transition into a continuous one as one increases 
the concentration of impurities. Albano 9 simulated the ZGB model on incipient percolation 
clusters (IPC's) with a fractal dimension of 1.90. He showed that both transitions, at y\ 
and y 2 are continuous, and that for an infinite lattice, in which y co is larger than 0.408, 
the reactions stop at finite times because the IPC's are poisoned by pure CO. Casties et 
al. 10 also performed a Monte Carlo simulation of the CO oxidation on probabilistic fractals. 
They observed a change in the character of the transition at y 2 from first order on regular 
lattices to second order on percolation clusters (for p larger than p c = 0.593, which is the 
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percolation threshold on the square lattice). 

In this work we have performed mean-field ( site and pair approximations) calculations 
and Monte Carlo simulations for different values of the concentration of inert sites. The 
model studied here is a variant of the original ZGB model, where inert sites are randomly 
distributed over the lattice. Our approach is close related to that presented by Vigil and 
Willmore 11 to study the effects of spatial correlations on the oscillatory behavior of a modified 
ZGB model, where defects are continually added and desorbed from the surface. In their 
studies, they considered the mean-field site and pair approximations, as well as Monte 
Carlo simulations. In the present work we have determined the phase diagram for different 
concentrations, and the spinodal and transition lines as a function of the concentration 
of inert sites. We have constructed hysteresis curves to find the critical concentration at 
which the first-order transition changes into a continuous one. This paper is organized as 
follows: in Sec. II we present the results obtained within the site approximation; in Sec. 
Ill we introduce the pair approximation equations and show the results obtained using this 
scheme; Sec. IV presents the results of simulations, and finally, in Sec. V we present our 
conclusions. 



We take a square lattice as our catalytic surface. A fraction rid of the sites is randomly 
distributed over the lattice representing the pinned inert sites. The remaining sites of the 
lattice can be vacant, or occupied by either O atoms or CO molecules. The ZGB model is 
described by the following steps: 



II. SITE APPROXIMATION 




(1) 



2 (g) + 2v^ 20(a), 



(2) 




(3) 
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where the labels g and a denote the gaseous phase and an adsorbed reactant on the surface, 
respectively, and v indicates a vacant site. Steps (1) and (2) indicate the adsorption of the 
species, whereas the third step is the proper reaction, between distinct species located at 
adjacent sites of the lattice. In the site approximation the time evolution equations of the 
concentrations are given by 

= -y™n v {l - (1 - n G ) 4 ) + 2(1 - y co K(l - n co ) 3 , (4) 
= Vco n v (l - n ) 4 - 2(1 - y co )n 2 v (l - (1 - n co ) 3 ), (5) 

where n Q , n co and n v represent, respectively, the coverages of O, CO and blank sites in the 
lattice. y co gives the arrival probability of a CO molecule. In addition, there is the following 
constraint among the concentrations 

n co + n G + n v + n d = 1. (6) 

The steady-state solutions of the above system of equations are given by n v = 0, that 
corresponds to a poisoned surface, and 

2(1 - Voo) 

Inserting Eq. (7) into Eq. (4) we obtain an expression for the steady-state values of the 
concentration n co : 



(n co + y + n d ) 4 + (l-n co ) 3 -l = 0. (8) 

We exhibit in Fig. 1 a typical diagram for the coverages of CO, O and vacant sites 
obtained for n d = 0.2. This diagram was obtained by integrating the equations of motion 
for the n co and n Q concentrations, starting from an initial condition in which the number of 
empty sites is n v = 1 — n d . The site approximation does not give any continuous transition 
for all values of the concentration of inert sites. This was already pointed out by Dickman 6 
for the ZGB model without inert sites. We observe in Fig. 1, that the limit of stability of 
the reactive phase is y s = 0.467510, which corresponds to the spinodal point. Therefore, a 



reactive steady-state is found for all values of y co < y s . For values of y co > y s , the system 
becomes poisoned, with a large amount of CO and a small concentration of O atoms. The 
presence of O atoms in the region y co > y s is due to the inert sites that can block some 
oxygen, and also to the simplicity of the site approximation, which does not forbid the 
formation of O — CO nearest-neighbor pairs in the lattice. The tolerance of these O — CO 
pairs also explains the absence of the continuous phase transition, which is observed in the 
simulations. Fig. 2 is a plot of the solutions n co of Eq.(8) versus the parameter Y for 
different values of concentration rid of inert sites. We obtain two solutions, which we call n> 
and nf , that join together at the spinodal point. For instance, for rid = the value we find 
is Y = 0.638986, which furnishes the value y s = 0.561013. We also note in Fig. 2 that, at 
the spinodal point, the concentration of n co molecules remains the same irrespective of the 
value we choose for rid- This special value is n co = 0.1660. Then, the net effect of adding 
rid is to shift the curves horizontally. In this site approximation, solutions are possible only 
for values of n d < 0.638986. This happens because above this value the solution would 
correspond to the non-physical value y co < 0. So, the meaning of the two solutions in 
Fig. 2 is the following: the branch nf Q represents the stable steady-state solutions whereas 
the n> branch gives the unstable solutions. These solutions were obtained after numerical 
integration of the equations of motion for n co and n Q , starting from the state described by 
n v — 1 — ri d - For the initial condition n v — Y and n co larger than ra> the system evolves to 
the poisoned state. The initial condition n v = Y and n co less than n> drives the system to 
the lower stable reactive solution nf . 

Fig. 2 also shows that, as we approach the spinodal point for any value of rid, the 
region of stability becomes narrower. Then, we expect that for some value of y co < y s a 
first-order transition occurs, that is, the concentration n co must increase abruptally from a 
small value ( reactive phase) to a large value (poisoned phase). Unfortunately, we cannot 
use here the usual thermodynamic considerations based on the minimization of a suitable 
thermodynamic potential. In order to find this first-order transition we adopt the same 
kinetic criteria employed by Dickman 6 , which was borrowed from the work of Ziff et al 5 . 
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The phase transition was determined by choosing an initial state where half of the lattice 
was empty and the other half was completely filled with CO. In this work we choose as our 
initial state, to solve the equations of motion for n co and n Q , the values n v = n co = \{l — n<i). 
It is clear that this choice is not the same as that considered by Ziff et al., because we cannot 
discriminate which sites are empty or not. The phase boundary is defined at the special 
value 2/2 where the solution of the equations of motion changes from the reactive to the 
poisoned state as we vary the value of y co for the same initial condition, as established 
above. For n d = we obtain the same value found by Dickman. We exhibit in Fig. 3 
the results obtained for the first-order transition and the spinodal points for rid = 0.2. The 
spinodal was obtained from the initial condition n v — 1 — rid, and the first-order transition 
from the condition n v = n co — |(1 — n d ). For this particular value of rid = 0.2, we have 
%)2 = 0.3999 and y s = 0.4821. We have considered all values of the concentration of inert 
sites, and Fig. 4 shows the values of y s (dashed line) and y 2 (full line) as a function of the 
concentration of inert sites. At the particular value rid = 0.55 the two lines merge. For 
values of rid > 0.55 the transition still remains of the first-order type, although the number 
of vacant sites that stay in the active state is very small. For instance, for n d = 0.60, at 
the transition point (y 2 = 0.0725), the number of vacant sites changes from 0.0388 in the 
active state, to 2 x 10~ 8 in the poisoned state. Throughout our analysis we considered a 
given state to be active if the number of vacant sites is larger than 10~ 6 . We also exhibit in 
Fig. 5 the number of vacant sites n v at the active state as a function of the number of inert 
sites rid, at the transition and spinodal points. We observe that for all values of rid < 0.55, 
the number of vacant sites at the spinodal point is always larger than that at the transition 
point. 



III. PAIR APPROXIMATION 

Let us consider the application of the pair approximation procedure to this ZGB model 
that includes inert sites. Here we introduce the pair probability P a p of a random nearest 
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neighbor pair of sites being occupied by species a and (3. We have the following types of 
species: v, d, c, and O, which represent, respectively, vacant, inert, carbon monoxide, and 
oxygen. As in the previous treatments 12,13 we need to consider only the changes that occur 
at a particular central pair in the lattice. In the table below we display the allowed and 
forbbiden (indicated by x) nearest-neighbor pairs in the present model. 





V 





c 


d 


V 


vv 


vo 


vc 


vd 


o 




oo 


X 


od 


c 




X 


cc 


cd 


d 








dd 



The next table also exhibits all the possible transitions among pairs. We obtain 14 
independent transitions, labelled by numbers in the range 1 — 14. In the table transitions 
indicated by x are prohibited. 



from > 

to I 


vv 


vo 


vc 


vd 


oo 


od 


cc 


cd 


VV 




4 
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X 


12 


X 


VO 


1 




X 


X 


10 


X 


X 


X 


VC 


2 


X 




X 


X 


X 


13 


X 


vd 


X 


X 


X 




X 


11 


X 


14 


oo 


3 


5 


X 


X 




X 


X 


X 


od 


X 


X 


X 


8 


X 




X 


X 


cc 


X 


X 


7 


X 


X 


X 




X 


cd 


X 


X 


X 


9 


X 


X 


X 





Then, we write the equations relating the probability of each element with the corre- 
sponding pair probabilities: 



8 



Py Pvv Pvo H~ Pyc H~ Pydi 

Pq Pod H~ Pvo H~ Poo ; . . 

(9) 

Pc Pcd ~\~ Pvc ~\~ Pec i 
Pd = Pdd + Pvd + Pod + Pcd- 

The pair probabilities also satisfy the constraint 

Pvv + Poo + Pec + Pdd + 2(P 1)0 + P vc + P^i + P d + Pcd) = 1. (10) 

Next, we need to write the time evolution equations for the pair probabilities. Examining 
the latter table we can construct the desired equations of evolution. We explicitly write the 
equations of motion for the pair probabilities P a p. 

dP 

til vo 



dt 
dP 

111 vc 

dt 

dPyd 

dt 
dP 

u 1 OO 

dt 
dPod 
dt 

dPcc 
dt 
dPcd 

dt 
dP 

dt 
dPdd 



tl + tlQ — £4 — £5, 

*2 + ^13 — ^6 — hi 
til + tlA — ^8 ~~ ^9) 
^3 + 2^5 — 2t 10 , 

ts-tn, (11) 

2tj — ti2 — 2ti3, 

t% — ii4, 

— £3 + #12 ~~ 2£i + 2^4 — 2t2 + 2^6, 

0. 



dt 

where t\ to t u are the transition rates. The factors of two arising in the equation of motion 
for P vv are due to the fact that the pair probabilities P^ and Pji are equal by symmetry. 
For instance, from the pair vv we can obtain, with the same probability, the different con- 
figurations ov and vo. In general, the expressions for the transition rates are lengthy, and 
we present these transition rates in the Appendix. 

In this pair approximation we cannot obtain analytical solutions as we have done in the 
site approximation. We solved the coupled set of eight nonlinear equations by the fourth- 



order Runge-Kutta method, searching for the stationary solutions. We considered the two 
different initial conditions as in the case of the site approximation. 

Let us first consider the evolution from the initial state where P v — 1 — P d , in which 
only the pairs vv, vd and dd are present in the lattice at t — 0. Figure 6 shows the diagram 
of the model for P d = 0.2. For < y co < 0.2299 the lattice poisons with oxygen. In the 
range 0.2299 < y co < 0.4821 there is an active region, and for y co > 0.4821 the lattice 
poisons with CO. When Pd = 0, we found the same figures obtained by Dickman in his 
pair approximation. For instance, the site and pair approximations give the same value 
for the spinodal point y s . However, when we consider some inert sites in the lattice, the 
spinodal point found in the site approximation is always smaller than that obtained within 
the pair approximation. For this particular value, P d = 0.2, the site approximation yields 
y s = 0.4675, whereas y s = 0.4821 is obtained by the pair approximation. The value of y co 
at the continuous transition, which now arises in this pair approximation, decreases slightly 
with increasing values of the concentration of inert sites. 

We also considered the solutions evolving from an initial condition where half of the free 
sites (P v = 1 — P d ) is filled with CO molecules and the other half left empty. In order to be 
close to the initial condition used in the simulation, we chose for the initial pair conditions 
Pec = P vv , P dv = Pdc and P vc = 0, which mimics a division of the lattice into two parts: 
on one side of the lattice we would have inert sites and CO molecules and, on the other 
side, vacant and inert sites. If P d = 0, we found for the transition between the active 
and CO-poisoned states the value y-i = 0.5240, which agrees with the value found in the 
simulations. Fig. 7 displays the concentration of CO molecules at the transition point for 
which P d = 0.2. In this pair approximation, the values of y s and yi are very close. We also 
show in Fig. 8 the concentration of vacant sites as a function of the concentration of inert 
sites, at the transition point, and also at the spinodal point. Both curves join at Pd = 0.50, 
and for P d > 0.60, we cannot observe any active state. As in the site approximation, an 
active state is defined only if P v > 10~ 6 . Then, the calculations performed within the pair 
approximation give results that are very similar to those obtained by the site approximation, 
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concerning the spinodal and transition points. 

In addition, it was observed that initial conditions do not affect the point in which 
the continuous phase transition occurs. In Fig. 9 we exhibit the phase diagram for this 
ZGB model with inert sites. The size of the reactive window decreases as we increase the 
concentration of inert sites. We have plotted the transition line for the first-order transition 
and for the spinodal line, which gives the limit of stability of the reactive phase. The line 
separating the active and O-poisoned phases is a continuous transition line. 

IV. SIMULATIONS 

We have performed Monte Carlo simulations in the ZGB model with inert sites in order 
to check the results we have obtained in the site and pair approximations. The simulations 
were carried for different values of the concentration of inert sites P d . For small values of 
P dl we considered square lattices of linear size L = 64, but for large values of P d we have 
used lattices of linear size up to L = 150. The first step in the simulation is to randomly 
distribute the selected fraction P d of inert sites in the lattice. All simulations then started 
with a fraction of empty sites equal to P v = 1 — P d . The CO molecules arrive at the 
surface with a probability y co and the 2 molecules with probability 1 — y co . The rules for 
adsorption and reaction of the species are exactly the same as in the original ZGB model 5 . 
Since adsorption of oxygen requires two nearest neighbor empty sites, the effect of the inert 
sites is to favour the adsorption of CO relatively to that of 2 molecules. In general, we 
have taken 10 3 Monte Carlo steps (MCs) to attain the stationary states, and 10 3 more to 
calculate the concentration averages at the stationary states. One MCs is equal to L x L 
trials of deposition of species, where L is the linear size of the lattice. To speed up the 
simulations we worked with a suitable list of empty sites. 

We exhibit in Fig. 10 the phase diagram of the model in the plane y co versus P d . It is 
similar to that obtained within the pair approximation. However, there is a fundamental 
difference between the transition line separating the active and CO-poisoned phases in both 
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approaches. In the pair approximation the transition line is always of the first-order type, 
whereas in the simulations there is a critical concentration above which the transition be- 
comes continuous. We have done detailed simulations to find the critical concentration at 
which the transition becomes continuous. We have found for the critical concentration of 
inert sites the value P% = 0.078. We arrived at this value by looking at the hysteresis loops 
in the curves of P co versus y co for different values of the concentration P d , as we can see in 
Fig. 11. We proceed as follows: in Fig. 11a we fixed the concentration of inert sites at the 
value 0.070 and the curve with circles, which is the proper transition curve, was obtained 
from an initial state where P v — 1 — Pd, that is, with a lattice almost empty. The curve 
with squares was determined from an initial state in which the lattice was almost covered 
by CO. We have taken a fraction of only 5% of randomly empty sites over the lattice at the 
starting time. Then, we clearly observe the hysteresis loop at the concentration P d , which 
implies that the transition is of first-order. On the other hand, Fig. lib, where the fraction 
of inert sites is P d = 0.080, does not exhibit the hysteresis loop and the transition is clearly 
a continuous one. The critical value of P% = 0.078 was obtained analysing the behavior of 
these curves in the range 0.070 < P d < 0.080. As we have pointed out in the Introduction, 
Hovi et al. 7 had already observed the change in the nature of this transition as a function 
of the concentration. The phase boundary separating the active and the O-poisoned phases 
in Fig. 10 is continuous for all values of P d . We have checked this fact by observing that no 
hysteresis loop was found for any value of P d . The width of the active phase decreases with 
increasing values of P d . For values of P d > 0.45 the lattice is poisoned (absence of empty 
sites) with different amounts of CO and O species. Due to finite size effects, this value is 
larger than the value 0.408 found by Albano 9 in the limit of very large IPC's. 

We have also noted that the production rate of CO2 molecules attains its maximum 
value exactly at the first-order transition, for values of P d < P d . If P d > P d the maximum 
production rate of CO2 molecules is located inside of the reactive window. This is seen in 
Fig. 12, where the circles indicate the points where the production rate of CO2 is maximum. 
In the site and pair approximations this maximum occurs always at the phase boundary, 
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irrespective of the value of Pa- Fig. 13 shows the production rate R of C0 2 molecules as a 
function of P d . As expected, the role of inert sites is also of blocking the reactions over the 
catalyst. The maximum production rate occurs at a surface free of impurities. 

V. CONCLUSIONS 

We have studied the effects of a random distribution of inert sites on the phase diagram 
of the ZGB model. We determined the time evolution equations for the concentrations 
of the different species over the catalytic surface within an effective field theory, at the 
level of site and pair approximations, and also performed Monte Carlo simulations on the 
model. We obtained the coverages of the species as function of the deposition rate of CO 
and of the concentration of inert sites. In the site and pair approximations we found the 
transition line and the limit of stability of the reactive phase. In the site approximation, the 
continuous transition between the O-poisoned and reactive states is absent for any values of 
the concentration of inert sites. The width of the reactive window exhibits the same behavior, 
as a function of concentration of inert sites, in both pair approximation and Monte Carlo 
simulations. However, the transition between the reactive and CO-poisoned phase is always 
of first-order in the site and pair approximations, whereas Monte Carlo simulations give a 
critical point where the transition changes nature. For values of the concentration of inert 
sites less than the critical value, the transition is first-order and above this value, it changes 
to a continuous one. The determination of this critical concentration was possible through 
the analysis of the hysteresis curves for different values of the concentration of inert sites. 
The production rate of CO2 molecules is maximum at the first-order transition, in both site 
and pair approximations. This is the case in the simulations, but the transition is of the 
first-order type. When the concentration of inert sites is greater than the critical value, 
the maximum production rate of CO2 molecules moves towards the reactive window. The 
overall effect of inert sites is to reduce the production of C0 2 molecules. 
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APPENDIX: TRANSITION RATES IN THE PAIR APPROXIMATION 

We present the transition rates in the pair approximation, which we used in Section III 
to solve the time evolution equations for the pair probabilities. The transition rates are t\ 
to ti4, which are given by 
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FIGURES 

FIG. 1. Coverages of CO (full line), O (dashed line) and empty sites (dotted line), as a function 
of the deposition rate of CO. The concentration of inert sites is rid = 0.2. Coverages obtained 
through the site approximation. 

FIG. 2. Stability curves in the site approximation for different values of concentration of inert 
sites. The upper (n> ) and lower (n< ) branches give the concentration of CO, respectively, at the 
unstable and stable states. The open circles indicate the position of the spinodal points. From 
the outer to the inner curves the concentration of inert sites is 0.0, 0.1, 0.2, 0.3, 0.4, 0.5 and 0.6. 

V — Sfco 

2(l-y co ) • 

FIG. 3. Plots of the concentration of CO (n co ) versus y co to locate the spinodal and transition 
points for the particular value rid = 0.2. Spinodal (dashed line), first-order transition (full line). 

FIG. 4. Behavior of the spinodal (dashed ) and first-order transition (full) lines as a function 
of concentration of inert sites in the site approximation. 

FIG. 5. Concentration of empty sites at the spinodal (dashed line) and at the first-order tran- 
sition (full line) as a function of the concentration of inert sites. 

FIG. 6. The same legend as in Figure 1, but coverages obtained within pair approximation. 

FIG. 7. The same legend as in Fig. 3. Calculations performed in the pair approximation. P co 
has the same meaning that n co . 

FIG. 8. Concentration of vacant sites at the spinodal (dashed line) and at the first-order tran- 
sition (full line) as a function of the concentration of inert sites, in the pair approximation. 
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FIG. 9. Phase diagram of the ZGB model with inert sites. Pd gives the concentration of inert 
sites, y co is the probability that the CO molecule hits the surface. We have an active phase and a 
O-poisoned and CO-poisoned phases. The dashed line represents the spinodal points, the full line, 
the first-order transition points, and the dotted line gives the continuous transition points. Phase 
diagram obtained in the pair approximation. 

FIG. 10. Phase diagram in the plane Pd versus y co , obtained through Monte Carlo simula- 
tions showing the reactive window. The squares give the transition between the reactive and 
CO-poisoned phases. The circles represent the continuous transitions between the active and 
O-poisoned phases. The lines are a guide to the eyes. 

FIG. 11. Hysteresis curves near the critical value of the concentration of inert sites, (a) 
P d = 0.070, (b) P d = 0.080. 

FIG. 12. Plot of Pd versus y co showing the points where the production rate of CO2 is maximum 
(circles), and the production (squares) at the transition between the active and CO-poisoned phases. 
The lines serve as a guide to the eyes. 

FIG. 13. Maximum production rate of CO2 molecules (R) as a function of the concentration of 
inert sites Pd- 
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